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ABSTRACT 

Fast brightness variations are a unique tool to prohe the innermost regions of active galactic nuclei (AGN). These variations are called 
microvariahility or intra-night variability, and this phenomenon has heen monitored in samples of hlazars and unobscured AGNs. 
Detecting optical microvariations in targets hidden hy the obscuring torus is a challenging task because the region responsible for 
the variations is hidden from our sight. However, there have been reports of fast variations in obscured Seyfert galaxies in X-rays, 
which rises the question whether microvariations can also be detected in obscured AGNs in the optical regime. Because the expected 
variations are very small and can easily be lost within the noise, the analysis requires a statistical approach. We report the use of a 
one-way analysis of variance, ANOVA, with which we searched for microvariahility. ANOVA was successfully employed in previous 
studies of unobscured AGNs. As a result, we found microvariable events during three observing blocks: in two we observed the same 
object (Mrk 477), and in another, J0759-I-5050. The results on Mrk477 confirm previous findings. However, since Mrk477 is quite a 
peculiar target with hidden broad-line regions, we cannot rule out the possibility that we have serendipitously chosen a target prone 
to variations. 

Key words, galaxies: active - galaxies: individual (Mrk 477, SDSS J075-t5050, SDSS J1430+1339), galaxies: quasars: general 


1. Introduction 


The different spectral features of AGNs can be explained by the 
different orientation of the object along the observer’s line of 
sight. Type 2 AGNs are oriented edge-on, and therefore the view 
of the central engine is blocked by the obscuring torus. Type 2 
AGNs can furthermore be divided according to their luminosity. 
While Seyfert 2 galaxies are low-luminosity type 2 AGNs, type 
2 quasars are the high-luminosity counterparts. Determining the 
luminosity of obscured AGNs is no trivial task because the dust 


obscuring torus effectively blocks our view. Therefore Reyes et 


al. ( 200 8 1 proposed a proxy, where obscured AGNs with L[oin] S 
10® are classified as type 2 quasi-stellar objects (QSOs). 


Microvariahility (or intra-night variability) is a phenomenon 
of variations with a low amplitude (hundredths of a magni¬ 
tude) over a short timescale (minutes to hours). The first re¬ 
ports of microvariations were made shortly after the discovery 
of quasars ([Matthews & Sandage||196^, but the reported mi¬ 


crovariations were considered unreliable because of technical 
difficulties. The first reliable detections of microvariations were 
reported by [Miller et al.j ( [1989| l in hlazars, which are extremely 
powerful radio sources. 

Microvariahility must originate in a region of some light 
minutes in size to agree with causality arguments. We can as¬ 
sume that this region is located near the center associated with 
the most violent physical phenomena. Microvariahility may arise 
either from instabilities in the jet or in the accretion disk, or both 
( Ramirez et al. [20091 2010| l, and therefore provides a unique op¬ 
portunity to study the unresolved inner regions. These regions 
are hidden in type 2 QSOs, and the contribution of the host 
galaxy to the total flux cannot be neglected, thus these objects 
are unlikely to show microvariations. 

This paper is organized as follows; in Sect. 2, we describe the 
sample, observations, and data treatment. Section 3 is dedicated 
to the analysis of the data and the brief description of the statis¬ 
tical test used. In Sect. 4 we summarize the results and discuss 
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Table 1. Properties of the observed targets from the catalog of |Reyes ^ 
[aTlpOOS} . Right ascension and declination are in J2000. Brightness in 
V filter was transformed from SDSS filters. The values correspond to 
the catalog of quasars and active nuclei: 13th edition (|Veron-Cetty &| 
|Veron|2010| ). 


Target 

R.A. 

dec 

^mag 

z 

J0759H-5050 

07:59:40 

-h50;50:24 

16.59 

0.06 

J1430H-1339 

14:30:29 

-h13:39:12 

17.64 

0.09 

Mrk 477 

14:40:38 

-h53:30:19 

15.03 

0.04 


them along with possible explanations in Sect. 5. At the end, we 
provide short conclusions. 


Table 2. Colors of the targets, reference stars, and comparison stars for 
every observed field. The values are adopted from the SDSS data release 
12 . 


Field 

star ID 

9 

g-r 

J0759H-5050 


16.09 

0.31 


starl 

15.58 

0.44 


star2 

13.65 

1.18 

J1430H-1339 


15.80 

0.37 


starl 

16.68 

1.1 


star2 

15.01 

1.21 

Mrk477 


14.84 

0.14 


starl 

15.01 

0.60 


star2 

13.23 

-1.05 


2. Sample selection, observations, and data 
reduction 

2.1. Sample 

Our aim was to obtain observations of a homogeneous set of 
obscured type 2 quasars. We have chosen type 2 quasars over 
Seyfert 2 galaxies because the contamination of the host galaxy 
is lower than in the Seyfert 2 galaxies as a result of the higher 
intrinsic brightness of the nuclei in type 2 QSOs. The sample 
was selected on the basis of brightness. It was chosen from the 
catalog provided by Reyes et^ (2008i, which is based on the 
Sloan Digital Sky Survey Data Release 6 ( |Adelman-McCarth^ 
|et al.|[2008] l. This catalog consists of 887 type 2 quasars se¬ 
lected based on [OIII] emission line luminosity. It is the most 
complete collection of type 2 QSOs to date. The limiting red- 
shift of the catalog is z = 0.83. The catalog covers 80% of the 
SDSS DR6 spectroscopic database. Our targets were chosen to 
be bright (g < 17 mag) to achieve accurate photometric mea¬ 
surements with the telescope. There are four objects fulfilling 
our requirements in the catalog. All of them are observable from 
the northern hemisphere. We were able to observe three of them. 
The fourth target met with telescope-pointing problems. The list 
of the observed targets and their properties can be found in Table 
[T] All the chosen targets are at z < 0.1. Thanks to this, we can 
compare similar properties between the sample objects because 
the same phenomenon falls into similar broadband filter. 


2.2. Observations and data reduction 

Our sample was observed between March 28 and March 31, 
2011, with the SITe 4 detector mounted on the 1.5 m telescope 
at San Pedro Martir observatory in Baja California, Mexico. We 
employed the same observing strategy as in previous works (de 
Diego et al. 1998[ Ramirez et al. 2009| l. All targets were ob¬ 
served in the Johnson V filter. Exposure time for each target was 
set to 60 seconds, providing a singal-to-noise ratio (S/N) > 100. 
This corresponds to a precision of 0.01 mag. The detector was 
binned (2x2) to allow fast readout with low noise. We observed 
all the targets with less than 1.3 airmass, and all of our observing 
nights were dark. We checked the colors of the stars that we have 
used for the analysis. [Carini et ar] ( |1992| l and |Stalin et al.|p()04| l 
proposed that the airmass effect is negligible as long as the color 
difference is g - r < 1.5. We provide the color information for 
our targets and stars in the Table 

Bias subtraction and flat-fielding were performed using the 
standard tasks from the IRAF data reduction softwar^D Cos- 


Image Reduction and Analysis Facility (http://iraf.noao.edu/l. 
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were removed using L.A. Cosmic]^ by 


van Dokkum 


mic rays 

(2001 . Aperture photometry was performed using SExtractor 
2.8.^ (Bertin & Arnouts 19961. We used fixed apertures for both 
the target and the stars in the field. We made sure to enclose all 
stars and host galaxies detectable in our images. This is impor¬ 
tant because the contribution can be variable for bright galaxies 
as a result of intra-night fluctuations in the atmospheric seeing, 
as pointed out by Carini et al. ( 1992|l. According to Cellone et 


al. P000|l, an aperture that would include the whole host galaxy 


can account for the seeing variations. For stars, we used aper¬ 
tures of 25 px (corresponding to 6.32 arcsec) diameter, while for 
the AGNs we used a 40 px (10.12 arcsec) diameter aperture. This 
is roughly four times the full width at half maximum (FWHM) 
of the objects. Hence we assume that with the apertures that we 
have chosen, the contribution of the host galaxy is constant over 
time and therefore is not causing any additional variations. 


3. Analysis 


We took advantage of the presence of other stars in the field to 
perform differential photometry. We used one star as a reference 
and another one for comparison to ensure the stability condi¬ 
tions. The left columns of Fig. 1 show the stability of the field 
stars during the observations. 

According to de Diego ( |2010| l, there are two statistical tests 
suitable to search for microvariability, namely an F-test and a 
one-way analysis of variance; ANOVA. ANOVA surpasses F- 
test in power in general, but an enhanced F-test ( de Diego|2014| l 
can provide more statistical power than ANOVA. Unfortunately, 
the enhanced F-test takes advantage of the presence of many ob¬ 
jects in the same field as the target to provide calibrating stars. 
Figure 13 of |de Diego| ( |2014| l shows that ANOVA is the most 
suitable test for our needs. 

The observing strategy was specifically chosen with the aim 
to use ANOVA. For each object, we divided the data from every 
night into groups of five data points, following the strategy of 
|de Diego et al.| ( [T998[ |. We computed the variance in each of the 
groups and the variance among the groups. We largely benefited 
from the fact that ANOVA does not rely on the error estimation 
from the photometric software. The errors are computed from 
the standard deviation within each group of five observations. 
Light curves with error bars estimated from ANOVA for every 
first observing night of each target are provided in Fig. [T] 


^ Laplacian cosmic ray identification/ 

(http://www.astro.yale.edu/dokkum/lacosmic/1. 

^ Source Extractor (http://www.astromatic.net/software/sextractorl. 
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Table 3. Results of microvariability in type 2 quasars. The second column describes the difference we used to create the differential light curve, 
where mgso stands for the target, stands for the reference star, and stands for the comparison star used to control the stability of the 
stellar light curve. N stands for the number of groups used as an input for ANOVA. In the hfth column, we provide the probability of obtaining the 
result if the null hypothesis is conhrmed. The level of signihcance is or = 0.001). 


Target 

DLC 

Date 

N 

Probability 

Microvariability (a - 0.001) 

J0759H-5050 

niQso - 

-rriref 

Mar 28 2011 

30 

0.0975 

N 


^comp 

- niref 

Mar 28 2011 

30 

0.132 

N 


rriQso - 

-mref 

Mar 30 2011 

20 

0.073 

N 


^^comp 

- ttlpef 

Mar 30 2011 

20 

0.146 

N 


mgso - 

- mref 

Mar 31 2011 

30 

5.5E-04 

Y 


^^comp 

Itlpef 

Mar 31 2011 

30 

0.327 

N 

J1430h-1339‘' 

mQso - 

-mref 

Mar 28 2011 

25 

0.0105 



^comp 

- mref 

Mar 28 2011 

25 

0.027 



mQso - 

-mref 

Mar 30 2011 

12 

0.391 

N 


^^comp 

- ITlfgj 

Mar 30 2011 

12 

0.436 

N 

Mrk 477 

mgso - 

- mref 

Mar 29 2011 

30 

lE-04 

Y 


^^cotnp 

Itlpef 

Mar 29 2011 

30 

0.114 

N 


rriQso - 

-mref 

Mar 31 2011 

20 

4.95E-05 

Y 


^comp 

- mref 

Mar 31 2011 

20 

0.102 

N 


Target J1430+1339 has the same probabilities for stars and target, all of them low. The lack of stars in the field does not allow 

conclusions based on this observation. 



Fig. 1. Differential light curve with errors computed as standard devi¬ 
ations within the sets of five observations for one of the targets and a 
comparison star. A constant is subtracted from the star to fit it into one 
plot. 

4. Results 

We have detected microvariability in two out of three targets that 
we observed in the run. It has been shown that there is a strong 
evidence, at the a = 0.001 level of significance, of variability in 
the light curves of threee observations of type 2 quasars. Below, 
we discuss the results for individual targets. We present light 
curves for the microvariable events in Fig. For the microvari¬ 
able events, we also provide the microvariability amplitudes as 
defined by |Heidt & WagneT| ( |1996| l, 

Amp — ^{Af^ax ~ 2cr^. (1) 

J0759-I-5050 was observed during three nights. The observa¬ 
tions from March 30, 2011, was disturbed by clouds, so it is not 
strictly continuous. 

Microvariability was only detected during the observations 
carried out on March 31. The event was observed during the pe¬ 
riod between 3:21 and 5:42 UT. The upper level of the amplitude 
is 0.04 mag, which is well above our estimated precision set by 
the high S/N. The variability amplitude is 4.6%. Because the mi¬ 
crovariations were only observed during one of the nights of the 


campaign, we conclude that we cannot confirm microvariability 
in this source. 

J1430-I-1339 was observed on March 28 and March 30, 
2011. During the two observing nights, both the target and 
the comparison star exhibited very contradictory results from 
ANOVA. We were unable to find a comparison star that was not 
variable at a = 0.001, although the p - values were likewise too 
low to confirm variability. Unfortunately, the field around this 
particular target limited the choice of comparison stars we could 
have used for ANOVA. In conclusion, there is not enough evi¬ 
dence of microvariability in this source. 

Mrk 477 was observed on March 29 and March 31, 2011. 
The gap between the observations on March 29 was caused by 
a technical maintenance pause. Mrk 477 is the only target that 
showed microvariable behavior in both observing nights. The 
upper limit for variations is 0.027 between 5:55 UT and 7:02 UT 
and 0.023 mag between 8:46 UT and 10:02 UT (for the first and 
second observing night, respectively). Computed amplitudes of 
microvariability are 2.7%and 2.3% for the fisrt and second ob¬ 
serving night, respectively. Based on our results, we conclude 
that Mrk 477 exhibits microvariations. 

5. Discussion 

We have detected microvariable events with a confidence level 
of a = 0.001 during three observing blocks. The positive de¬ 
tections were made for J0759-H5050 (one observing block out of 
three) and Mrk 477 (microvariability detected in two out of two 
observing blocks). J0759H-5050 was only variable during one of 
the observing blocks, we were unable to confirm the variability 
behavior in the other observing nights. The amplitudes of the 
variations as computed using the approach by |Heidt & WagnCT 
\\99b\ are comparable to the other quantitative measurements of 
microvariability ( [Romero et al.|1999] l. It is highly unlikely that 
our measurements are affected by systematics because possible 
effects would also affect the comparison stars, which we have de¬ 
tected as not variable. Because we observed during dark nights, 
we can rule out any variation caused by significant change in the 
background level. We also controlled the airmass of the obser- 
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Fig. 2. Differential light curves for Mrk 477 (panels a — d) and J0759+5050 (panels e — /) for the nights with detected variability. Panels a, c, and 
e represent niQso — rtiref, while panels b, d, and / represent nicomp - iriref- The reference star light curves do not show any features, as expected. 


vations. Since the airmass does not rise above 1.3 airmasses, it 
is highly unlikely that any of the variations would be based on 
the color change of the source. Given that the target detected as 
variable is in the vicinity of other sources in the image, we con¬ 
structed check graphs with different aperture sizes to validate 
that there is no spurious light from other sources propagating 
into our measurements. The apertures are sufficiently large to be 
unaffected by a change of the seeing. 

Mrk 477, the brightest and closest of our sources, is the only 
one that was widely observed before our study. It was already re¬ 
ported to show optical microvariability by |Jang| ( (2001| l, although 
the report lacks details and comments on robustness of the result 
and the statistics used. |Tran| ( |1995b| l reported detection of hid¬ 
den broad-line regions (HBLR) as a part of spectropolarimetric 
observations. |Tran] ( | 1995b) l found that the degree of polarization 
is P = 1.2% and the rotation of the position angle is © = 85°. 
The presence of HBLR might provide a clue about the mecha¬ 
nism behind the variability of the source. Until now, almost all 
the sources that show optical microvariability are type 1 AGNs, 
meaning that they have a visible broad-line region, and we can 
assume that part of the variations originates there. Our detection 
of variability would therefore be consistent with the detection 
of the HBLR reported by |Tran| ( |l995b[ ). The mechanism that al¬ 
lows us to detect the variability remains unclear, although there 
are possible explanations. 

It is possible that the observed light originates close to the 
central engine and is then reflected on the corona of charged par- 
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tides above the dusty torus. Polarimetric observations would be 
needed to confirm such a scenario. We expect that microvaria¬ 
tions in obscured sources show small amplitudes, therefore this 
scenario would agree with our results. Nevertheless, we cannot 
confirm this statement as we would need to observe the source 
for an extended amount of time in polarized light. 

Mrk 477 is also the only target out of our sample that was 
detected in X-rays. It was reported as a Compton-thick source 
( |Gilli et al.|20I0| l. Recently, there have been reports of spectral 
variability in X-ray detected obscured sources (most notably the 
Seyfert 1.8 galaxy NGC 1365; [RTsaliti et aL]|2009[ [Brenneman| 
|et al.|201'3) on short timescales. The variations are caused by al¬ 
ternate covering and uncovering of the X-ray source by clouds 
in the broad-line region. Even though there are archival data of 
our source from the XMM-Newton mission, the monitoring fre¬ 
quency is insufficient to perform a study of short-term X-ray 
variations. Therefore we cannot infer whether similar behavior 
occurs in Mrk 477. Nevertheless, we cannot rule out this possi¬ 
bility. 

Another possible explanation lies within the structure in the 
dusty torus, which could provide observational windows to ob¬ 
serve the inner parts of the otherwise obscured AGNs. The orig¬ 
inal unified scenario introduced by Antonucci (19931 includes a 
geometrically thick obscuring rotating region, or in other words, 
a torus. It is worth noting that |Krolik & Begelman ( 1988) 1 sug¬ 
gested that such a structure would be hard to maintain and intro¬ 
duced a clumpy structure into the model of the torus. The model 
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introduced by [Nenkova et al.| ( |20081 l counts no more than five 
or ten clouds, which are unsuitable to explain high-frequency 
variations. Considering orbital motion for the individual clouds 
in the dusty torus, variability is expected whenever these clouds, 
temporarily, do not completely hide the central engine. This may 
be the case for a few Seyfert galaxies that have shown spectral 
variations, changing from type 2 to type 1 AGNs and vice versa. 
However, these variations occurred on timescales of years, and 
our data do not allow us to observe such a behavior. 

6. Conclusions 

Microvariation in type 2 quasars is a phenomenon previously 
unaccounted for. From our observing campaign, we reported a 
detection of one variable (Mrk 477) and one unconfirmed vari¬ 
able type 2 quasar (J0759H-5050) that were detected as variable 
during one observing block shorter than two hours. Both results 
were obtained using a statistical method, namely a one-way anal¬ 
ysis of variance, ANOVA, which has previously proved to be a 
powerful test for the detection of microvariations in AGN light 
curves. Given the observation history of the variable object Mrk 
477, we discussed the possible causes of the variations. Our ob¬ 
serving campaign was limited to one filter in the optical regime, 
therefore we cannot confirm any of our hypotheses by obser¬ 
vational evidence. Because of the possible importance of mi¬ 
crovariability to the study of the physics of the type 2 quasars 
and their obscured regions, it is worth investigating whether this 
is a ubiquitous phenomenon in this class of AGNs. 
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